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Abstract. The Hanle effect in strong spectral lines is the physical mechanism that 
should be increasingly exploited for quantitative explorations of the magnetism of the 
quiet solar chromospheric plasma. Here we show, by means of multilevel radiative 
transfer calculations and new spectropolarimetric observations, that the amplitude and 
shape of the scattering polarization profiles of the Ha line is very sensitive to the 
strength and structure of the chromospheric magnetic field. The interpretation of the 
observations in terms of one-dimensional radiative transfer modeling suggests that there 
is an abrupt magnetization in the upper chromosphere of the quiet Sun. 



1. Introduction 

Spectroscopy of the hydrogen Ha Une at 6562.8 A has played a crucial role to visualize 
that the "quiet" solar chromosphere is a fibrilar dominated- magnetism medium (e.g., 
the review by Rutten 2007), but quantitative information on the magnetic field vector 
can only be obtained through spectropolarimetry. 

Unfortunately, the spectral line polarization produced by the Zeeman effect is of 
little practical interest here because the magnetic splitting of the cr and tt components of 
the line transitions are practically negligible compared with the very significant Doppler 
width of the Ha line, so that in the quiet Sun the Zeeman effect only produces observ- 
able circular polarization signals of very low amplitude (with V//~0.01% when using 
today's telescopes), which nonetheless are practic ally insensitive to th e physical condi- 



tions of the true chromospher e (see the review by Trujillo Bueno| l2010l) . 



Fortunately, as shown by lStepan & Trujillo Bueno (2010), the amplitude and shape 
of the Ha linear polarization profiles produced by scattering processes in the solar at- 
mosphere are very sensitive (through the Hanle effect) to the strength and structure of 
the magnetic field in the upper chromosphere of the quiet Sun. Their conclusion was 
based on radiative transf er modeling of th e peculiar line core asymmetry (LCA) of the 
Qjl profile observed by iGandorfeij fcOOd) in a quiet region close to the solar limb, an 
asymmetry which is however not seen in the observed I{A) profile. The main aim of 
this paper is to clarify further why we believe that the measurement and interpretation 
of the QII and U jl profiles produced by scattering processes and the Hanle effect in 
Ha might lead to an important breakthrough in our empirical understanding of chro- 
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Figure 1. Left panel: Grotrian diagram of neutral hydrogen showing Ha (red 
lines), Lya (blue lines), and Ly/3 (violet lines). The four red solid lines correspond 
to the transitions which make a significant contribution to the solar Ho- polarization. 
Right panel: Zeeman splittings of the FS levels of « = 3 calculated in the incomplete 
Paschen-Back effect regime. The red arrows indicate the Bfj values of the individual 
levels. 



mospheric magnetism. To this end, here we show both additional radiative transfer 
calculations and an example of the spectropolarimetric observations we are obtaining 
with the Zurich Imaging Polarimeter (ZIMPOL) attached to the Gregory Coude Tele- 
scope (GCT) of the Istituto Ricerche Solari Locarno (IRSOL). Interestingly, these new 
observations show non-zero signals in both Q/I and U/I. 



2. Basic information about the Hanle effect in Ha 

The physical origin of the so-called scattering line polarization is nothing but the selec- 
tive emission and absorption processes produced by the presence of atomic alignment 
in the line's levels, which results from radiative transitions produced by anisotropic 
radiation (e.g., Trujillo B ueno 2001). 

The Ha line is produced by the 3£j-2fj' transitions of neutral hydrogen (see 
the left panel of Fig. [T]l. The Zeeman splitting of the upper levels of Ha is shown in 
the right panel of the same figure. The selection rule Al = ±1 for radiative transi- 
tions inhibits radiative couplings between levels (n, I, j) and {n, l±\, /). Therefore, 
for B<100 G we can safely neglect quantum coherences between different j-levels. 
The thermal broadening of the Ha line is large due to the high chromospheric temper- 
atures and low atomic mass of hydrogen. A typical Doppler width is about 0.3 A. For 
magnetic strengths B<100 G the wavelength shifts of the Zeeman cr and n components 
do not produce any noticeable contribution to the Stokes Q and U profiles, which are 
dominated by the presence of atomic level alignment and its modification by the Hanle 
effect. Since the separation of the fine structure (FS) components of the Ha line is about 
0.1 A, the ensuing transitions do overlap. 

Via the Hanle effect each of the FS levels is sensitive to a diflFerent strength of 
the magnetic field according to the expression Bh = 1.137 x 10"^/fiife^> where Bh 
(in gauss) is the so-called critical Hanle field of the level under consideration, tufe its 
lifetime in seconds, and g its Lande factor. The linear polarization of the scattered 
radiation is mostly sensitive to the Hanle effect of magnetic fields with intensities of the 
order of Bh- If the magnetic field is weaker than about 0.1 Bh the polarization state of 
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Figure 2. Left panel: Line source functions in a semi-infinite isothermal atmo- 
sphere with «f, = Hp = 4 X 10'" cm and T = 10"* K. Right panel: The corresponding 
fractional anisotropy, Ji/J^, of the line radiation. 



the scattered radiation is not modified. If B > \00 Bh, the polarization state is sensitive 
only to the field direction (i.e., the so-called Hanle saturation regime). The critical fields 
of the upper levels of Ha range from approximately 6 to about 15 G. Note also that only 
the FS levels with angular momentum j > 1/2 may carry atomic alignment and thus 
contribute to the scattering polarization of the Ha line. 



3. Radiative Transfer Calculations 

Our radiative transfer calculations of the Ha scattering polarization profiles and their 
modification due to the action of the Hanle efi"ect are carried out within the framework 
of the density-matrix theory of spectral U ne polarization, which uses the approxim ation 
of complete frequency redistribution (see lLandi Degrinnocenti & Landolfi 12004 ). The 



numerical solution of this type of problem requires finding, at each spatial grid point of 
the model atmosphere under consideration, and for each j-level of the chosen atomic 
model, the self-consistent values of the multipolar components of the density matrix. 
We achieve this through the application of efficient iterative meth ods and formal solvers 
of the Stokes-vector transfer equations (see lTrujillo Buenoll2003L for a detailed review). 



3.1. Ha in an isothermal unmagnetized atmosphere 

We assume a 3-«-level model atom consisting of 9 FS levels and consider a semi- 
infinite, isothermal model atmosphere with temperature T = 10^ K and constant elec- 
tron and proton densities (i.e., Uf. = Up - 4 x 10^^cm~^). T he inelastic coUisional rates 
of the transitions are taken from iPrzybilla & Butlei^ ( 2004 ) and th e depolarizing colli 



sional rates are calculated in the semi-classical approximation of iSahal-Brechot et al. 
(11996). We have calculated the self-consistent solution of the ensuing non-LTE scat- 
tering polarization problem. Fig. |2] shows the line source functions and the radiation 
anisotropics of the Lyor, Ly^S, and Ha lines throughout the atmosphere. 

In order to calculate the emergent linear- polarization, one has to take into account 
the atomic polarization of the FS levels. In the left panel of Fig. [3] we show the frac- 
tional alignment of the individual levels as a function of height in the atmosphere. At 
the formation depths of Ha, the strong resonant Lya line is optically thick and the Lya 
radiation is virtually isotropic there (see the right panel of Fig. Consequently, the 
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Figure 3. Fractional level alignment (left panel) and the emergent fractional linear 
polarization profile for fi - Q.l (right panel) without including any magnetic field in 
the isothermal atmosphere. The arrows indicate the atmospheric heights where the 
line optical depth is unity for line-of-sights (LOS) with yU = 1 and /i = 0. 1 . 




Figure 4. Same as Fig. [3] but for a magnetic field with a fixed inclination Ob - 65° 
and with a random azimuth changing its intensity with height in the atmosphere. 



polarization of the 2/73/2 level is negligible and, therefore, selective absorption pro- 
cesses in Ha can be neglected. The Ha scattering polarization is fully determined by 
polarization of the upper levels, especially by the 3dj ones. Thanks to this fact, the 
emergent Q/I profile can be expressed as the sum o f the contributions of the indivi dual 
FS transitions (see the right panel of Fig. 12 cf . IStepan & Trujillo Buenol llOlOl for 
details). 



3.2. Ha in an isothermal magnetized atmosphere 

For simplicity, here we consider the case of a magnetic field with a fixed inclina- 
tion Ob with respect to the vertical and with a random azimuth, xb, uniformly dis- 
tributed in the interval [0, 27r]. Such a magnetic field does not break the cylindrical 
symmetry o f the atmosphere and makes the numerical solution relatively simple (see 



Stepan & Trujillo Buena.201 1,. for details). If the magnetic field intensity is uniform 



throughout the atmosphere, the fractional level alignments Pq/Pq are reduced as the 
field i ntensity increases and the amp litude of the Q/I line profile is reduced accord- 
ingly (IStepan & Truiillo Buenoll201lh . 
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Figure 5. The emergent Q/I profiles at /i - 0.1 for different azimuths of the 
horizontal magnetic field (see the different line styles and the;^'^ values in degrees). 
Left panel: Magnetic field intensity is decreasing with height around t - I. Right 
panel: Magnetic field intensity is increasing with height around r = 1 . 



If there is a spatial gradient in the magnetic field strength, then the shape of the 
resulting 2/7 profile may be more complicated than the typical Gaussian bell-Iike shape 
of Fig. [3] Such an effect is due to the different magnetic sensitivities of the (blended) 
Ha line components and to the fact that different parts of the line profile originate at 
slightly different depths in the model atmosphere. A clear demonstration of this effect 
can be achieved for the case in which the magnetic field intensity increases with height 
around the height corresponding to unit line optical depth. An academic example of 
such a magnetic field variation is shown in the left panel of Fig. \4\ The line center 
photons originate in the upper layers of the atmosphere which are strongly depolarized 
by the magnetic field there, and the line-center amplitude thus decreases. Given that the 
different components of the line are contributing differently to the resulting Q/I profile 
and since they are sensitive to different field strengths, a line-center asymmetry (LCA) 
may appear (see the right panel of Fig. 



3.3. Ha in a semi-empirical model with magnetic field gradients 

In Fig. [5] we show examples of synthetic Ha Q/I profiles obtained by as suming two 



types of magnetic field gradients in the semi-empirical model atmosphere of Fontenla et al 
(11993, hereafter, the FAL-C model). While the model with magnetic strength decreas- 
ing with height produces quite symmetric Q/I profiles for all the field azimuths con- 
sidered, the model with magnetic strength increasing with height leads to a flattening 
of the Q/I profile and to a LCA which can be significantly stronger than when assum- 
ing a random-azimuth field, especially for xb ~ 0° and 180°. The corresponding U/I 
profiles do not show any noteworthy asymmetry and they are thus not plotted here. 
Calculations for the case of a r andom-azimuth field can be found in Figs. 2 and 3 of 
IStepan & Truiillo Buenol ( l2010l) . 



3.4. Ha in hydrodynamical models of cbromospheric dynamics 

The solar chromosphere is a dynamical system whose time evolution is far too fast to 
be resolved by spectropolarimetric observations. Consequently, only the time average 
of the emergent line profiles can be measured with today's telescopes. Fortunately, 
radiative transfer synthesis of the emergent scattering polarization in the available ID 
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Figure 6. Left panel: Time evolution of the temperature (solid line) and of 
the proton densi t y (da sh-dotted line) in one of hydrodynamical simulations of 
ICarlsson & SteinI (119971) . Right panel: The corresponding amplitude of the Ha frac- 
tional linear polarization emerging at yU = 0. 1 (solid line). The horizontal dashed line 
shows the amplitude of the time averaged Q/I profile. 



hydrod ynamical models of chromospheric internetwork dynamics by ICarlsson & Stein 
19971) are a suitable step prior to performing similar calculations in 3D models of the 



quiet chromosphere. 

W e haye used the hydrodynarn ical simulation referred to as the strongly dynamic 
case by lAsensio Ramos et al.l (l2003l) and we have solved the Ha scattering polarization 
problem in each ID model corresponding to each time step of the hydrodynamical sim- 
ulation. The results can be found in Fig. |6] The left panel shows the time variation of the 
kinetic temperature and proton density (which is the main source of collisional depolar- 
ization of Ha) at the height where the line optical depth for a line of sight with = 0. 1 is 
unity. One can identify in Fig.[6]the shock waves going through the atmospheric region. 
These are followed by a slow recombination of hydrogen. In the right panel of the same 
figure, the amplitude of the emergent Q/I profile at // = 0.1 is plotted as a function of 
time. Note that after the wave passes through the line formation region (which abruptly 
decreases the polarization signal), the emergent linear polarization slowly increases un- 
til the next shock arrives. These time intervals of increasing amplitudes correspond to 
the compression phases of the simulated chromospheric plasma. The Q/I profiles with 
the largest amplitudes are thus red-shifted in the spectrum while the blue-shifted phases 
correspond to the smallest signals following immediately the shock wave passing. In- 
terestingly, the average 2/7 profile, whose amplit ude is very similar to the one ob tained 
in the non-magnetic FAL-C model (see Fig. 2 of (Stepan & Trujillo Buenoll2010h . does 
not show any noteworthy asymmetry at all. 



4. Spectropolarimetric observations and modeling 

The suggestion by lStepan & Trujillo Buenol (|2010|) of the existence of an abrupt change 
in the degree of magnetization of the upper chromosphere of the quiet Sun was based on 
ID radiative transfer modeling of the peculiar line core asymmetry of the Q/I profile 
observed by Gandorfer (|2000|) . using semi-empirical models of the solar atmosphere 
such as the FAL-C model. Given that the real solar chromosphere is highly inhomoge- 
neous and dynamic any such ID static model can only be considered as a poor repre- 
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Figure 7. An example of our recent spectropolarimetric observations of the Ho- 
Stokes profiles at about 8 arcsec inside the solar limb (dashed lines) obtained with 
ZIMPOL at IRSOL and the theoretical fit (solid lines) we have obtained by assum- 
ing an abrupt magnetization and an enhanced microturbulent velocity in the upper 
chromosphere of the FAL-C model. 



sentation of the cbromospheric thermal and density stratification. Therefore, we cannot 
exclude the possibility of an alternative explanation. Nevertheless, it is intereting to 
note that the temporally averaged Q/I profile calculated without magnetic fiekis in the 
hydrodynamical simulations of chromospheric dynamics of ICarlsson & SteinI (119971) 
does not show any noteworthy LCA (see Section IX4l ). 

There are two i nter-related avenues to confirm o r refute the above-mentioned ten- 
tative conclusion by [Stepan & Trujillo Buenol (|201C1|) . On the one hand, it is impor- 
tant to develop realistic three-dimensional MHD simulations of the quiet solar chro- 
mosphere, and to u se them to compu te the emergent Q/I profile for comparison with 
the observations of iGandorferl ( 2000l) . On the other hand, it is crucial to carry out full 
Stokes observations of the He line at various on-disk positions. We are presently work- 
ing on both fronts, but here we focus only on showing an interesting e xample of one of 
our re cent spectropolarimetri c observations carr ied out with ZIMPOL jGandorfer et al. 



20041) at the GCT of IRSOL (Bianda etaL||2009'), which in addition to showing a similar 



LCA in the observed Q/I profile it shows a clear non-zero U/I signal (see the dotted 
lines of Fig. |7]l. The observations were obtained near the heliographic North pole at 
jU = 0.13 ± 0.01 on June 13th 2009. The spectrograph slit width was set to 0.5 arcsec 
and the total exposure time was 800 seconds. The data redu c tion i ncluded the correc- 
tion for instrumental polarization described by iRamelli et all (l2006l) . The Q/I and U/I 
images show interesting structures and variability along the spatial direciton. The pro- 
files reported here have been obtained integrating 16 arcsec along the spectrograph slit. 
The solid lines of Fig. |7] show a good fit to these new observations, which we have ob- 
tained by using the FAL-C model assuming a n abrupt magnetization abo ve a height of 
1900 km (i.e., similar to that of Fig. 3A in Stepan & Trujillo Buenol (1201 0)) and a fixed 
magnetic field inclination {9b~65°) and azimuth (^5*135°). It is however important to 
mention that although the core of Q/I is well fitted when using the height-dependent 
microturbulent velocity of the FAL-C model, we had to artificially increase it up to 
about 20km s~' in order to be able to achieve a better fit of the profile wings. We inter- 
pret this as an indication of broadening by unresolved motions and/or structures in the 
real solar chromosphere. 
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5. Conclusions and Perspectives 

The scattering polarization of the Ha line and its sensitivity (via the Hanle effect) to 
magnetic fields between and 50 G is a promising tool for diagnostics of chromo- 
spheric magnetic fields. In this paper, we have described the peculiar character of the 
fine structure components of the Ha line and of the formation of its scattering po- 
larization profiles in optica lly thick model atmosph eres, including the time-dependent 



mospr i 

hydrodynamical models of ICarlsson & Steiiil (11997 ). Moreover, we have presented an 



interesting example of the full Stokes vector observations we are carrying out, which 
shows a conspicuous LCA in Q/I and a clear U/I signal. 

For future progress, new observations of the Ha polarization with higher spatial 
resolution within and outside coronal holes are urgently needed. The modeling of such 
observations with the help of the new generation of chromospheric models that are 
currently under development will facilitate new significant advances in the exploration 
of chromospheric magnetism via spectropolarimetry. 
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